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Abstract 

In this paper we present an efficient numerical approach based on the Renor- 
mahzation Group method for the computation of self-similar dynamics. The latter 
arise, for instance, as the long-time asymptotic behavior of solutions to nonlin- 
ear parabolic partial differential equations. We illustrate the approach with the 
verification of a conjecture about the long-time behavior of solutions to a certain 
class of nonlinear diffusion equations with periodic coefficients. This conjecture 
is based on a mixed argument involving ideas from homogenization theory and 
the Renormalization Group method. Our numerical approach provides a detailed 
picture of the asymptotics, including the determination of the effective or renor- 
malized diffusion coefficient. 
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1 Introduction 



The time evolution of many non-equilibrium physical systems is well described by 
nonlinear partial differential equations (PDEs). It is frequently observed that under 
suitable, but fairly general, conditions the time evolution of such systems becomes 
asymptotically self-similar. Within the PDE framework, this assertion amounts to the 
observation that the PDEs in question have solutions that behave asymptotically as 



It turns out that often such a self-similar behavior possesses a great deal of universal- 
ity: the scaling exponents a and (5 and most aspects of the profile function (/)(■) are 
independent of initial conditions (within a suitable class) or, more strikingly, even of 
the form of the equations. 

Universal behavior is a central issue in the study of critical phenomena in equilibrium 
statistical mechanics and quantum field theory. Using the Renormalization Group 
(RG) approach [H], HSl [Hj, physicists predict critical exponents and determine the 
universality class of a variety of models. In the early 90 's, Goldenfeld, Oono and 
collaborators (see [12] and references therein) developed a perturbative renormalization 
group method for PDEs and used it to compute perturbatively the similarity exponents 
in the difficult and most frequently encountered case of self-similarity of the second kind 
(in Barenblatt's classification [Ij). Later, Bricmont, Kupiainen and Lin [SIE] introduced 
a non-perturbative renormalization group approach. A numerical renormalization group 
algorithm was developed at the same time by Chen and Goldenfeld [8]. 

In this paper we combine the numerical approach of Chen and Goldenfeld with the 
constructive approach of Bricmont et al. to develop a computationally efficient RG 
algorithm for the calculation of self-similar dynamics. This algorithm is then used 
to produce a detailed picture of the asymptotic dynamics of certain nonlinear diffusion 
equations with periodic coefficients; in particular, this algorithm allows the investigation 
of which aspects of these dynamics are universal and which are not. A similar point of 
view was adopted by some of the authors [1] to study a modification of Barenblatt's 
equation. Isaia [T^ also used the same approach to study other classes of problems 
including some with logarithmic corrections to the time decay. 

Our purpose in writing this paper is threefold. Firstly, we conjecture the long-time 
asymptotic behavior of a class of PDEs. Secondly, following the strategy of [U [13] , we 
present a systematic numerical approach, based on the RG ideas, to study self-similar 
dynamics. Finally, we point out a connection between the renormalization group and 
homogenization theory [2]. Homogenization theory has been applied with success to 
critical lattice field theories [16], where a RG approach should work. Here, in addition 
to adopting a reverse point of view and showing that the RG approach succeeds in 
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problems typically analyzed with homogenization theory, we further point out that the 
RG approach gives more general results; see below. 

Numerical procedures based on rescaling, and thus similar in spirit to the RG approach 
presented here, have previously been developed [1313] and used to study solutions which 
blow up in finite time (as, for example, in the case of certain react ion- diffusion equations 
and the cubic nonlinear Schrodinger equation). Such procedures exploit the known 
self-similar structure of the solutions under study to define the appropriate rescalings. 
Recently, versatile and more efficient versions of such procedures were developed [T71 [TU] 
and employed for the computation of solutions which blow up at multiple points. Yet, 
the RG procedure presented here is unique in exploiting fixed points. The central 
feature of this procedure is that the rescalings are determined dynamically, not by any 
a priori knowledge of the scale invariance of the solutions. As a result, the procedure can 
be used to investigate self-similar asymptotics of both the first and second (anomalous) 
kind. On the other hand the current implementation of the procedure that we present 
here is not appropriate for studying blow-up problems. 

This paper proceeds as follows: in Section [2] we define the class of initial value problems 
we study and state the conjecture we make about the long time behavior of their 
solutions. In Section [3] we present a non-rigorous argument to motivate this conjecture. 
In Section H] we explain how we have implemented our modification of the numerical 
renormalization group scheme of Chen and Goldenfeld [8] to the oscillating equation ([1]). 
In Section O we present numerical results that corroborate our conjecture. We also 
present studies to validate our numerical procedure, including the computation of the 
critical exponent a for solutions of Barenblatt's equation to illustrate the versatility 
of our procedure in computing anomalous behavior. Finally, in Section [6] we present 
numerical evidence supporting a similar conjecture for "relevant" perturbations (see the 
Remark in the next section), make some comments and discuss further developments. 

2 Asymptotics of oscillating equations 

Consider the following initial value problem (IVP): 



m(x,1) = /(x), X G M. 

Here /i and A are real parameters; [i is such that (1 + iigix)) > for all a; G M; g{x) is a 
smooth periodic function with zero mean; f{x) is a smooth and rapidly decaying initial 
condition; F{u, v, w) is an analytic function of u, v and w around u = v = w = 0: 




X e M, t>l 







{i>a, j>b, k>c} 
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Let T be the period of the function g and let H{g,fi) denote the harmonic mean of 
l + fig{x): 



H{g,fi) = 
We conjecture the following. 



-, -1 



T Jo l + 



dx 



(2) 



Conjecture 2.1 Let u{x, t) be the solution of the IVP (Qp. If the lowest order exponents 
of the power series expansion of F satisfy the inequality a + 26 + 3c — 3>0 and if the 
initial condition f is small enough, then there is a constant A, which usually depends 
on f , g, fi, X and F, such that 

A ^2 

t°'u{Vix,t) —>■ — e~4^ as t ^ oo, (3) 

v47r(T 

with a = 1/2 and a = H{g, /x). 



As we shall see, this conjecture arises from formal renormalization group arguments 
applied to ([1]). These arguments are based on the methods developed by Goldenfeld, 
Oono and collaborators [T2] and by Bricmont, Kupiainen and Lin |6l[7]. Implementing 
a modification of the RG-based numerical algorithm of Chen and Goldenfeld j8j, we 
verify the validity of the conjecture. 

Remark: Let dp = a + 26 + 3c — 3. Bricmont et al. [7] classify the perturbations F (see 
Eq- dl])) as irrelevant if dp > 0, marginal if dp = {F = uux, Uxx or u^) and relevant if 
dp < Q {F = Ux or m"). Conjecture 12. II is stated for irrelevant perturbations but similar 
conjectures hold for the marginal and relevant cases, although care is needed in order 
to avoid blow up at finite time. In particular, if \F = — irrelevant perturbations, 
1 < a < 3, then, instead of ([3]) we expect that 

t~u{\/tx,t) ^ fa{—^) as t — > oo, (4) 
where fa is a function such that 



u{x, t) = 




is a scale invariant solution to Eq. ([T]) with = and XF = — m", see [SI El Ej- In 
Section [6] we present numerical results supporting this conjecture. Due to logarithmic 
corrections, our numerical approach can not handle marginal perturbations efficiently. 
We point out that logarithmic corrections have been obtained using another version of 
the numerical RG in |13j . 

We emphasize the following aspects of ([3]): 



4 



1. the exponent a is universal: a = 1/2 independently of the initial condition /, the 
periodic function g, the parameters and A and the perturbation F{u,Ux,Uxx) 
provided that a + 26 + 3c — 3>0; 

2. the effective diffusion coefficient cr depends on fi and g: a = H{g,^); 

3. in general, the prefactor A depends on the initial condition /, on the perturbations 
F{u,Ux,Uxx) and g{x) and on the parameters fi and A. 

Conjecture 12.11 can be understood as a statement about the universality class of Gaus- 
sian fixed points of a renormalization group transformation. It furnishes a classification 
of the nonlinear terms that can be added to the linear oscillating equation without 
changing the long-time asymptotics of the associated solutions. The work of Duro and 
Zuazua [H] can also be interpreted within this framework. They have considered the 
equation 

ut - div(a(x)VM) = d- V{\u\'''\) in x (0, oo), 

where d G and a{x) is a smooth, symmetric and periodic matrix. Using homogeniza- 
tion theory, they have proved, among other results, the asymptotic result ([3]) whenever 
q > 1 + 1/N. Although we do not prove it, our conjecture indicates that Duro and 
Zuazua's results should hold for other non-hnear equations. 



3 Renormalization group analysis 

In this section we present the formal RG arguments which led to the aforementioned 
conjecture. In the sequel, we assume that the perturbation F{u,Ux,Uxx) in Eq. ([T]) is 
of the form u"'u\u'^xx- Later, it will become clear that this assumption constitutes no 
loss of generality since the higher powers of F are "irrelevant" in the RG sense. The 
renormalization group approach that we employ in this paper is simply the integration 
of the equation followed by a rescaling. To explain this idea, we need some preliminary 
notions. 

Let M be a real- valued function of (a;, t) G M x M+. For a fixed L > 1 and sequences of 
positive scaling exponents, {an}'^=i and {[3n}'^=i, define a sequence {u„}5^o rescaled 
functions inductively hy Uq = u and, for n > 1, 

Un{x, t) = (L^"x, Lt) = {L^"+-+^'x, L'^t) . (5) 

If the original function m is a global solution to Eq. ([1]), then a direct calculation reveals 
that Un satisfies the renormalized initial value problem: 

dtUn = Xn [1 + lig{uJnX)]dlUn + An U^^OxUnf {dlunf , (6) 
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Un{x, 1) = fn{x). 

Here 

_ rn-2(/3„+--+/3i) _ r/3n + -+/3i 

An — -'^ , l-i^n — -t^ ; 

and 

^ ^ j^^n-(6+2c)(/3„+--+/3i)+{l-a-fe-c)(a„+--+ai) ^-j?^ 

and the initial data /„ is 

fn{x) = L""+-+"^M L") . (8) 

The renormahzation group transformation of Bricmont et ah, which we denote as Rl, 
is now introduced. This transformation acts on the space of initial conditions and is 
defined having the linear diffusion equation 

ut = u^x (9) 

in mind. As explained in [SI Ej, the study of the long time asymptotics of solutions to 
([T]) with /i = is equivalent to studying the fixed points, and their basins of attraction 
(i.e. universality classes), of the transformation Rl. The RG transformation is designed 
so that its fixed points are the similarity solutions to ([9]). Below we adopt the same 
strategy and use Rl to study the long time asymptotics of solutions to the oscillating 
equation ([T]). 

To define Rl, take L > 1 and set ai = (3i = 1/2 for all i>l. With this choice, Xn = 
ujn = Xn = AL"[3-("+2b+3c)]/2 ^j^g jyp p becomes 

= [1 + fig{L-/'xMun + XL^l^-i^+'^+'^^y' <(5.n„)^(a^u„)^ (10) 

Unix, 1) = fnix). 

If Un{x, t) is the solution to the IVP (ITUI) for t G [1, L], Rl is defined as 

{RLfn){x) = L^'^Un {L^'^X,L) = /„+i(x). (11) 

Rl depends on n but we do not write explicitly this dependence to simplify the notation. 
Notice that Rl has the semigroup property: 

Rlf{x) = RL^fix) = L^'^u {L^/^x, L") , n = 1, 2, . . . . 

This fact allows us to investigate the limit ([3]), with t = L", by iterating Rl, and this 
is how we proceed. We start our analysis by reviewing the results of Bricmont et al. 
[HI [7] when n = 0^ X: 
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1. 



let TZ denote the linearized RG transformation (take /i = A = in equation ([T])). 
One can check that any multiple of the Gaussian distribution 



e 4 



'A-n 

is a fixed point of IZ. is called the Gaussian fixed point; 

2. any smooth function g decaying sufficiently fast and with 5f(0) = 0, where g{k) 
denotes the Fourier transform of g{x), contracts to zero under the action of TZ: 

C 

il^5'l|oo < -^Iklloo < llS'lloo, 

if L > C; 

3. any function / can be written as its projection onto 0=,, plus a remainder g with 
^(0) = 0: / = /(O)0, + ^7; 

4. write BJlf{x) as BJlf{x) = An(l)*{x) + gn{x), where gn{0) = 0. Then, if the 
exponents a, b and c are such that [(a + 26 + 3c) — 3] > 0, An —>■ A and gn{x) —>■ 
as n OO. 



In other words, Bricmont et al. proved that the long-time behavior of the solution 
u{x, t) to the IVP (lH) with /i = is given by the solution of 

u{x, 1) = A(f)^:{x). 

Now consider the case 7^ = A. Using homogenization theory [2J, one can conclude 
that the limit ([3]) exists. This is done as follows. After n RG iterations of equation ([1]) 
with A = we obtain: 

dtU=[l + fig{L^/^x)]dlu (12) 

Un{x, 1) = fn{x). 

Take e,^^ = L"/^. Observing that e„ ^ as n — > 00, we see that studying the asymp- 
totics of the solution to IVP (fT2l) as n — 00 is equivalent to studying the asymptotics 
as e„ — s> 0, a situation which can be dealt with by homogenization theory which says 
that the long time behavior is governed by a multiple of the solution to 

dtu = au^x (13) 

m(x, 1) = 0<^(x), 
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where a = H{g,fi) is the harmonic mean defined by equation ([2]) and 

^^{x) = -j==. 

V47r cr 

For the general case /i 7^ 7^ A, the RG iteration of the I VP ([1]) is represented by the 
IVP f lTU]) . As — > cxD, the second term on the right hand side of Eq. flTUl) is expected to 
be driven to zero if [(a + 26+3c) — 3] > 0. Writing R1f{x) as BJlf{x) = y4„0o.(x) + (y'„(a;), 
with ^"„(0) = 0, we also expect that An ^ A and that (y'„(x) — ^ as n — > 00. Therefore, 
for large values of n, we expect the solution to the IVP ([1]) to behave as the solution 
to the IVP ([13]), with A = limA„ and a = H{g,fi). 



4 The numerical procedure 

The RG-based numerical procedure for the integration of the IVP ([1]) is now introduced. 
The renormalization procedure constructs the sequence of renormalized IVP ([6]) (viz. 
the sequences {an}, {Pn} and {fn}) satisfied by the solution of ([T]) in the following 
steps. 

Start with /o = /, the initial condition of the IVP ([1]). For n = 0, 1, 2, . . . 

1. Evolve the initial data fn forward, from the initial time t = 1 to t = L, with 
L> 1 arbitrary but fixed, using Eq. ([6]), to obtain the solution Un{x,L). 

2. Compute the scaling exponents a„+i, Pn+i as 

MMl _ _^(o^ . _i .... 

- Un{0,L)- Un{0,Ly ^"+1-2- 
(This choice will be discussed shortly.) 

3. Define /„+i(x) = L''"+^Un{L^''+^x, L). 

Notice that if /o = /, then and the definition of /„ imply that 

fnix) = (L'^-x, L) = L""+-+°1m (L^^+ ' +z^ix, L") , (15) 

where u{-, L") is the solution of the initial value problem ([T]) at time t = L". 



A simple consideration of ([T5|l uncovers the relationship between the limiting behavior 
of {/n} as n — > 00 and the long time asymptotics of m. To wit, rewrite (|T5|) as 



where 

_ JJian-(an^ hQl) Q _ JJll3n~{l3n-\ h/Jl) 

Thus, if the hmits An ^ A, ^ B, an ^ <y, ^ P and /„ — as n — >■ oo are 
attained, we might expect that 

ina^finp^ L") = ^ J „ ( Bn , ,^ A(f){Bx) U OO. 

This hmit, in turn, estabhshes the self-similarity of the long time asymptotics of u. 

It is the numerical algorithm consisting of steps [T] — TSl above supplemented by the 
calculation of the prefactors {An} and {-Bn} that we use to verify the Conjecture 12.11 
about the asymptotic behavior of solutions to ([1]). We remark that this procedure 
yields a very detailed picture of the asymptotics. In particular, it seems ideal to study 
which aspects of the dynamics are universal and which depend on parameters and initial 
conditions. A number of comments are in order. 

The rationale behind the idea we use to compute {an} in step[2]is the reputed self-similar 
asymptotic dynamics we want to compute: in the self-similar regime, the solution u at 
a; = 0, behaves as 

M(0,t) ~ const t"". 

Thus, given the relation between Un in the time interval t E [^^L] and m in t G [L", 
(c.f. ([5])), we expect that, as n ^ oo, a„ approaches a. Whereas the determination of 
the exponents a„ can always be performed as in step [2] above, the determination of the 
exponents j3n is problem dependent. It usually involves a scaling relation between the 
exponents a„ and /5n so that certain (a priori chosen) parts of the differential operator 
remain invariant under the rescaling of step [31 To illustrate this point, consider the 
initial value problem ([T]) with F{u,Ux,Uxx) = u"'u\u%^. If /?„ is chosen to be 1/2, then 
the linear operator 

ut = u^x (16) 

remains invariant under the scaling defined in step [31 In this way, the dynamics associ- 
ated with f[T6l) can be explored. On the other hand, if Pn is chosen to satisfy the scaling 
relation (c.f. equation ([7])) 

1 - (6 + 2c)pn + {l-a-b- c)an = 0, (17) 
then it is the nonlinear operator 

Ut = Xu'^ulul, (18) 



that remains invariant. In this case, it is the dynamics of fll8l) that can be explored. In 
the RG language, the choice /5„ = 1/2 focuses the attention on the dynamics of f[T6l) and, 
as such, is suitable for the investigation of its irrelevant and marginal perturbations. 
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This is done in Section [5] for the Equation ([T]). When the perturbation is relevant, then 
a new scahng relation may be needed. For instance, with the choice f|T7j) . the nonlinear 
operator ut = u°'uJ'Uxx'^ remains invariant (under the associated scaling). Therefore this 
is the appropriate scaling (which then defines an appropriate RG map) if we want to 
study irrelevant and marginal perturbations to the nonlinear operator Ut = u'^uJ'Uxx^. 
In Section [H] we study a case when both the Laplacian and the non-linear term become 
marginal under an appropriate scaling (see the remark on Section [2] and the limiting 
behavior given by (jlj)). 

Finally, the spatial scaling in step [3] can be realized in two different ways: by rescaling 
the mesh size Ax without changing the discrete sites j = 0, 1, . . ., so that after one 
iteration the new mesh size is (Ax)i = L~^^Ax and the new mesh points are located 
at X = jE'^^^Ax; this is the approach adopted in Ref. or by rescaling the discrete 
sites while keeping the mesh size fixed, so that after one iteration the new discrete sites 
L~^j are located at x = L~^^jAx. In this approach, the values of the solution u at the 
fixed mesh points x = j'Ax have to be interpolated from the data given at the (new) 
discrete sites in each iteration. 

In our investigation of the asymptotics of problem ([T]) we have used the first approach. 
The successive spatial scalings performed in step 3 change the frequency of the periodic 
diffusion coefficient. It is actually this mechanism that leads to a renormalization of the 
diffusion coefficient in the long time regime. Thus, to resolve properly the effects of the 
oscillating diffusion, the mesh size has to be shrunk at the same rate as the frequency 
is increased. 



5 Numerical results 

For the discretization of the governing PDE, any appropriate scheme can be employed. 
We chose a simple explicit finite difference scheme that combines Euler's method for 
the time discretization with the standard three-point formula for the discretization of 
the Laplacian operator and centered differences for the first order spatial derivatives. 
Given the stability constraints, the resulting scheme is second order accurate. 

To eliminate the need of numerical boundary conditions in our computations of purely 
initial value problems with compactly supported initial data, additional grid points with 
associted zero-valued data were added at the boundaries of the computational domain 
before each timestep was performed. In all simulations 27 grid points were used for the 
discretization of the initial data in —5 < x < 5 (initial computational domain). The 
time step was dynamically chosen to satisfy the stability condition (1 + max \ fig\)At < 
CAx^ (the grid spacing Ax shrinks as the renormalization procedure is iterated). The 
constant C was chosen to be 0.45. 
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Table 1: Description of the numerical simulations, with L = 1.021 and (3 = 0.5 



n 


M 


9 


X 


a 


b 


c 


/ 


1 


.1 


91 











/l 


2 


.1 


91 













3 


-.15 


9i 











/l 


4 


.1 


91 











/3 


5 


.1 


91 


.1 


4 








h 


6 


.1 


91 


.1 


2 


1 





h 


7 


.1 


92 


.3 


8 








./a 


8 


.8 


91 











/i 


9 


.1 


93 











/i 


10 


.8 


93 











/i 


11 
















/i 


12 
















/2 


13 
















/3 


14 


.1 


91 


.1 


1 


1 


1 


h 


15 


.6 


93 


.1 





1 


1 


./2 



Table [T] summarizes our numerical simulations: the first column is the simulation num- 
ber; the other columns specify the parameters /i and A, the exponents a, h and c and the 
functions / and g appearing in the IVP ([6]) (functions fi and Qi, i = 1, 2, 3 are specified 
in Figures d] and [2]) . In all simulations, the scaling factor L was chosen to be L = 1.021. 
Although, in principle, any value of L > 1 can be used, for a too large value of L, a 
large number of time steps per RG iteration would be required, and the system would 
shrink too quickly (due to the spatial rescaling); for a very small L, a large number of 
RG iterations would be necessary to produce good enough accuracy. 




Figure 1: Initial conditions /i (continuous line), /2 (dashed line), /s (dotted line). 

We now discuss some of the simulations. In order to validate the numerical procedure 
and to illustrate its accuracy, efficiency and the amount of detailed information it 
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Figure 2: Periodic functions: gi{x) = cos(x), g2{x) = [cos(x) + sin(2x) + cos(4x)]/2.72, 
g3{x) = 1 — 2|x/7rmod2 — 1| 

provides, the algorithm was apphed to the heat equation (/x = = A), an equation whose 
asymptotic behavior is well known (simulations number 11, 12, 13). FigureslHUU and[5] 
refer to this study, as explained below. Also, simulations for Barenblatt's equation 



where H{u) is the Heaviside function, were performed for different values of e. Figure [6] 
shows the anomalous exponent a computed by our simulations, compared with the 
first-order approximation [12] 



Figure 3: Validation simulation (heat equation): plot of x n, for simulations 11, 12 and 
13, verifying that a* = lima„ is universal with respect to initial conditions. 

Figure [3] depicts the convergence of the computed exponent an to the theoretical value 
a = 1/2 as the numerical RG map is iterated for several distinct initial conditions. 



ut = {1 + eH{-ut))u, 
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Fi gure 4: Validation simulation (lieat equation): plot of A.^ X n, for simulations 11, 12 
and 13, verifying that An — > A and that A is mass dependent. 



Figure 5: Validation simulation (heat equation): plot of — /„_i||/||/n|| x n, for simu- 
lations 11, 12 and 13, verifying that the profile function converges upon the iteration 
of the numerical RG map. 



Figure 6: Relation between a and e for Barenblatt's equation. The continuous curve is 
the first-order approximation for a(e). The circle points are the results of our simulations. 

This illustrates the universality of a with respect to sufficiently localized (e.g. with 
compact support) initial conditions. Figure H] shows the convergence of the prefactors 
An to the theoretical values A (dotted lines). The latter depend on the mass of the 
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initial condition. 

We also verified that the profile function is equal to the Gaussian distribution, in agree- 
ment with the analytical results. Finally, in Figure \5\ we plot the relative difference 
between successive profiles, ||/n — /n-i| |/| |/n| I, as a function of the number n of RG 
iterations. The convergence of this difference to zero can be used as a diagnosis of 
the self-similar behavior of the long-time asymptotics; it can also be used to check the 
convergence of the procedure and as a practical stopping criterion for the number of 
RG iterations. 

The remaining figures. Figures [TlfTOl illustrate various aspects of the asymptotics of 
the solutions to the I VP ([I]) and corroborate our Conjecture 12. 1[ In all simulations 
performed, a„ was observed to converge quickly to the number of RG it- 

erations increases. In order to verify the universality of this fact, we varied: initial 
conditions (Figure [7]); nonlinear perturbations (Figure [8]); values of /x and periodic 
functions (Figure [9]). Figure [TO] shows the convergence of the prefactors An and their 
nonuniversality. We emphasize that in all simulations, the convergence of the relative 
difference ||/„ — /n-i||/||/n|| (in both and L°° norms) to zero was observed, as the 
number of RG iterations increased, thus providing strong evidence of the self-similar 
nature of the asymptotics. 




Figure 7: Plot of a„ x n, for simulations 1, 2 and 4, verifying that a is universal with 
respect to the initial condition. 

Finally, the effect of the term fig on the diffusion coefficient is illustrated in Figures [11] 
and [121 These figures display the points (— log^*, — logcp), where (p denotes the profile 
function /„ computed after the last RG iteration and 0*(a;) = exp(— x^/4). It is clear 
that for each simulation these points lie on a straight line. This implies that the 
computed profiles are the "renormalized Gaussians" 4>a{x) = exp(— a;^/4cr), with a 
being given by the slope of the corresponding straight line. The theoretical values 
for a are also plotted (stars) in these figures and they are obtained from the explicit 
computation of the harmonic mean of 1 + fig. 
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Figure 8: Plot of Q!„ x n, for simulations 1, 5, 6, 14 and 15, verifying that a is universal 
with respect to the addition of nonlinear terms. 



Figure 9: Plot of a„ x n, for simulations 1, 3, 7 and 9, verifying that a is universal with 
respect to the periodic perturbation g{-) and the parameter fi. 



Figure 10: Plot of An x n, for simulations 1, 2 and 3, showing that An does converge as 
n ^ oo. 
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Figure 11: Plot of — log((/)*) x — log{(j)), for simulations 8, 10 and 15. The stars plotted 
are the theoretically computed values for a. 




Figure 12: Plot of — log(0*) x — log((/»), for simulations 1, 3 and 11. The stars are the 
explicitly calculated a. 

6 Concluding Remarks 

From the numerical simulations we have a solid basis that indicates the correctness 
of our Conjecture 12. 1[ Furthermore, using the numerical algorithm of Section HI we 
have some preliminary results for the dynamics for the following equation (which is in 
divergence form): 

ut = dx{[l + fxg{x)]da:u) + \F{u, u^, Uxx) X G M, t > 1. (19) 

Equations ([T]) and (fTUj) differ by a multiple of g\x)ux. Our results indicate that this 
term is irrelevant in the RG sense even though the term is relevant. When F = u'^u^ 
in equation f|T9l) . our conjecture has been proved [9]. 

The version of the numerical procedure that we described here is not appropriate for 
the analysis of the dynamics when marginal perturbations of the type uu^ are included. 
Marginal perturbations have been handled successfully using a modified version of the 
procedure in |13]. 
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To verify the strength and flexibihty of the method, we also performed a numerical study 
of Equation ([1]) with XF = —u"', 1 < a < 3. First suppose = in Eq. ([1]). In this 
case, the Laplacian and the nonlinear term \F are equally important asymptotically. 
The choice /3„ = | keeps the Laplacian invariant to scaling. As discussed in Section HI 
we expect that a;„ will converge to a = as — > cxo, so that the nonlinear term 
remains in the limiting equation describing the long-time asymptotic behavior. This 
result is rigorous and it was first proved by Brezis, Peletier and Terman [5] and also 
verified by Bricmont et al. [Gj [7] using the RG approach. When /i 7^ 0, we expect that 
the same scalings will hold. These scalings lead to the conjecture (jl]), which we have 
verified numerically. In Figure [T3] we plot the curve a = together with the results 
of our simulations. The renormalization of the profile function by a factor of ^/a^ in 
the argument has been verified: we have plotted the computed profile function against 
x/a/o", where the factor a was computed a priori as the harmonic mean of 1 + fig{x), 
and the plot coincides with the plot of the computed profile function with fj, = 0. 



Figure 13: Relation between a and a for XF(u) = — u", with 1 < a < 3. The continuous 
curve is obtained theoretically and the stars are the results of our simulations. 
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